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ABSTRACT 

Using deep near-infrared imaging of the Hubble Deep Field South with ISAAC on the Very Large 
Telescope we find 6 large disk-like galaxies at redshifts z — 1.4 — 3.0. The galaxies, selected in K s 
(2.2/xm), are regular and surprisingly large in the near-infrared (rest-frame optical), with face-on effective 
radii r e = 0'.'65 - 0'.'9 or 5.0 - 7.5 kpc in a ACDM cosmology, comparable to the Milky Way. 
The surface brightness profiles are consistent with an exponential law over 2 — 3 effective radii. The 
WFPC2 morphologies in Hubble Space Telescope imaging (rest-frame UV) are irregular and show complex 
aggregates of star-forming regions ~ 2" (~ 15 hj kpc) across, symmetrically distributed around the 
-RT s -band centers. The spectral energy distributions show clear breaks in the rest-frame optical. The 
breaks are strongest in the central regions of the galaxies, and can be identified as the age-sensitive 
Balmer/4000A break. The most straightforward interpretation is that these galaxies are large disk 
galaxies; deep NIR data are indispensable for this classification. The candidate disks constitute 50% 
of galaxies with Ly > 6 x 10 10 /i 70 2 L Q at z = 1.4 — 3.0. This discovery was not expected on the basis 
of previously studied samples. In particular, the Hubble Deep Field North is deficient in large galaxies 
with the morphologies and profiles we report here. 

Subject headings: galaxies: evolution — galaxies: high redshift — infrared: galaxies 
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1. INTRODUCTION 

Disk galaxies are believed to undergo a relatively simple 
formation process in which gas cools and contracts in dark 
matter halos to form rotationally supported disks with ex- 
ponential light profiles (Fall & Efstathiou 1980; Mo, Mao, 
& White 1998). A critical test of any theory of galaxy for- 
mation is to reproduce the observed properties and evolu- 
tion of galaxy disks. 

Previous optical spectroscopy and HST imaging have 
yielded a wealth of data on disk galaxies at z < 1. (e.g., 
Vogt et al. 1996, 1997; Lilly et al. 1998; Barden et al. 2003), 
although contradictory claims have been made regarding 
the implications for the size and luminosity evolution with 
redshift (see Lilly et al. 1998; Mao, Mo, & White 1998; 
Barden et al. 2003), and the importance of surface bright- 
ness selection effects (see Simard et al. 1999; Bouwens & 
Silk 2002). 

It is still unknown what the space density and proper- 
ties are of disk galaxies at substantially higher redshift. 
Many galaxies at z ~ 3 have been identified using the ef- 
ficient U-dropout technique (Steidel et al. 1996a, b). Most 
of these objects are compact with radii ^1 — 2 hj^ kpc, 



while some are large and irregular (Giavalisco, Steidel, & 
Macchetto 1996; Lowenthal et al. 1997). However, the U- 
drop selection requires high far-UV surface brightness due 
to active, spatially compact, and unobscured star forma- 
tion. As a result, large and UV-faint disk galaxies may 
have been overlooked and, additionally, the morphologies 
of LBGs could just reveal the unobscured star-forming re- 
gions rather than the more evolved underlying population 
which forms the disk. 

The most direct evidence for the existence of large disks 
at high redshift has come from observations in the NIR, 
which provide access to the rest-frame optical. Here the 
continuum light is more indicative of the distribution of 
stellar mass than in the UV and nebular lines are accessi- 
ble for kinematic measurements, van Dokkum & Stanford 
(2001) discuss a iGselected galaxy at z = 1.34 with a ro- 
tation velocity of ~ 290 km s _1 . Erb et al. (2003) detect 
~ 150 km s~ x rotation at ~ 6 hjQ kpc radii in the Ha 
emission line of galaxies at z ~ 2.3 and Moorwood et al. 
(2003) find > 100 kins" 1 rotation at ~ 6 h^ kpc from the 
center of a galaxy at z = 3.2, seen in the NIR spectrum of 
the [O III]A5007 A emission line. 
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Fig. 1. — Left panels: The WFPC2 U300, averaged Vme + I$i4 (rest-frame UV), and our averaged ISAAC H + K s images 
(rest-frame optical), scaled proportional to Fx with arbitrary normalization per galaxy. The rest-frame UV morphologies are 
complex and symmetric with respect to the center of the smoother optical distribution. Middle panels: SExtractor's 7si4-band 
segmentation map, the smoothed /814 images, and the K s -baxiA images after subtracting the scaled, smoothed /814-band. From the 
segmentation map it follows that detection in 7~8i4 would likely split up most of the sources. The central residuals in K s — I smoot h 
demonstrate that optical and NIR light are distributed differently and that all galaxies have a "red" nucleus. Right panels: The 
radial profiles in the if s -band. The abscissa and the ordinate are respectively the mean geometric radial distance and the surface 
brightness along elliptical isophotes. The arrows mark la confidence intervals for measurements with signal-to-noise less than 1. 
Overplotted are the best-fit exponential law (dashed), r 1 / 4 law (dotted), and point+exponential (dash-dot) for galaxy 494 and 611. 



The imaging data in these studies, however, are of lim- 
ited depth and resolution, making it difficult to determine 
morphological properties. In this Letter, we present an 
analysis of the rest-frame ultraviolet-to-optical morpholo- 
gies and spectral energy distributions (SEDs) of 6 large 
candidate disk galaxies at z ~ 1.4 — 3 using the deepest 
groundbased NIR dataset currently available (Labbe et al. 
2003). Throughout, we adopt a flat A-dominated cosmol- 
ogy (Q M = 0.3, A = 0.7, H = 70 h 70 kms^Mpc- 1 ). All 
magnitudes are expressed in the Johnson photometric sys- 
tem. 

2. OBSERVATIONS 

We obtained 102 hours of NIR J S ,H, and K s imag- 
ing in the HDF-S (2'5 x 2f5) under excellent seeing 

7 http://www.strw.leidenumv.nl/~fires 



(FWHM«0"46), using ISAAC (Moorwood 1997) on the 
VLT. The observations were taken as part of the Faint 
InfraRed Extragalactic Survey (FIRES; Franx et al. 
2000). We combined our data with existing deep optical 
HST/WFPC2 imaging (version 2; Casertano et al. 2000), 
in the E/30O) -B450, V606 and Jsi4 bands, and we assembled 
a if s -selected catalog of sources with SExtractor (Bertin 
& Arnouts 1996). Photometric rcdshifts and rest-frame 
luminosities were derived by fitting a linear combination 
of empirical galaxy spectra and stellar population models 
to the observed flux points (Rudnick et al. 2001, 2003a). 
The reduced images, photometric catalog, and redshifts 
are presented in Labbe et al. (2003) and are all available 
on-line at the FIRES homepage 7 . Furthermore, we ob- 
tained optical spectroscopy with FORS1 on the VLT for 



Large disk-like galaxies at high redshift 



3 



Table 1 

Properties of high redshift disk galaxies in the HDF-S 



Galaxy 




z 


M c B , rest 


^0,B,rest 




T 1/2.K 


r i/2,I 




302 


19.70 


1.439* 


-22.70 


19.70 


0.89 


0.70 


0.86 


0.46 


267 


19.98 


1.82 


-22.88 


19.92 


0.75 


0.74 


0.88 


0.37 


257 


20.25 


2.027* 


-23.08 


19.53 


0.74 


0.74 


0.84 


0.36 


657 


20.68 


2.793* 


-23.56 


19.33 


0.76 


0.70 


0.74 


0.18 


611 


20.53 


2.94 


-23.59 


18.51 


O.Q5 j 


0.52 


0.97 


0.27 


494 


21.14 


3.00 


-23.31 


18.84 


0.75 J 


0.56 


0.86 


0.47 



a Catalog identification numbers (see Labbe et al. ! 

b A's-band total magnitudes 

c Rest-frame absolute B-band magnitudes 

d Face-on rest-frame B-band central surface bright 

e Face-on best-fit effective radii (arcsec) 



^ K a half-light radii (arcsec) 

9 7~8i4 half-ligt radii, PSF-matched to K s (arcsec) 
h Ellipticity 

1 Spectroscopic redshifts 

3 Two-component models (point + exponential) 



some of the sources (Rudnick et al. 2003b). Additional 
redshifts were obtained from Vanzella et al. (2002). As 
discussed in Rudnick et al. (2001, 2003a), our photomet- 
ric redshifts yield good agreement with the spectroscopic 
redshifts, with rms \z spec - z phot \/(l + z spec ) w 0.05 for 

Zspec ^ 1.4. 

Large disk galaxies in the HDF-South were identified 
by fitting exponential profiles convolved with the Point 
Spread Function (PSF) to the if s -band images. Six ob- 
jects at z > 1.4 have effective radii r e > 3.6 h^Q kpc, three 
of which have spectroscopic redshifts. The mean redshift 
of the sample is 2.4. We will focus on these large galaxies 
in the remainder of the Letter. The structural properties 
of the full if s -selected sample will be discussed in Trujillo 
et al. (2003). 

3. REST-FRAME OPTICAL VERSUS UV MORPHOLOGY 

The large galaxies are shown in Figure 1. They have 
a regular morphology in the ISAAC if^-band (2.2//m), 
which probes rest-frame optical wavelengths between 5400 
and 9000 A. In contrast, the WFPC2 1^06 and 7 8 i 4 -band 
morphologies, which map the unobscured star-forming re- 
gions at rest-frame UV wavelengths between 1500 and 
3300 A, are irregular with several knots up to ~ 2" (~ 15 
hjQ kpc) apart, symmetrically distributed around the K s - 
band centers. In a few cases the observed optical light is 
spatially almost distinct from the NIR. 

As a result of the structure in the WFPC2 imaging, 4 
of the large objects have been split up into two sources 
by Casertano et al (2003). Fig. 1 shows the correspond- 
ing "segmentation" map by "SExtractor" which illustrates 
how the pixels in each image are allocated to different 
sources. The galaxies were not split up when the K s band 
image was used to detect objects. However, the broader 
PSF in the if s -band image can play a role: if we smooth 
the /8i4-band data to the same resolution, wc find that 
SExtractor only splits up 1 galaxy. 

Hence, the question remains whether these 4 objects 
split up in 7814 are superpositions or whether they are 
part of larger systems. We tested this by subtracting the 
PSF-matchcd 7 8 i4-band images from the i^-band images. 
The /si4-band images were scaled to the if s -band images 



to minimize the residuals. The residuals are shown in Fig. 
1. In all cases, we find strong positive residuals close to the 
centers of the objects as defined in the X s -band, whereas 
residuals at any of the 7si4-band peaks might be expected 
in case of a chance superposition. Furthermore, we per- 
formed photometric redshift analyses for subsections of the 
images and found no evidence for components at different 
redshifts. 

4. profile fits and sizes 

Next, we fitted simple models convolved with the PSF 
(FWHM«0'/46) to the two-dimensional surface bright- 
ness distributions in the i^-band. The images are well- 
described by a simple exponential law over 2 — 3 effective 
radii (galaxy 302, 267, 257 and 657) or by a point source 
plus exponential (galaxy 611 and 494), where the point 
source presumably represents the light emitted by a com- 
pact bulge contributing about 40% of the light. We also 
derived intensity profiles by ellipse fitting. As can be seen 
in Figure 1, most galaxies are well described by an ex- 
ponential. The central surface brightnesses and effective 
radii, enclosing half of the flux of the model profile, are 
corrected to face-on and shown in Ta ble 1. The central 
surface brightness is multiplied with VI — e, as an inter- 
mediate case between optically thin and optically thick, 
and corrected for cosmological dimming. 

The effective radii (semi-major axes) are surprisingly 
large, r e = 0'.'65-0'.'9 (5.0-7.5 kpc in a ACDM cosmol- 
ogy), comparable to the Milky Way and much larger than 
typical sizes of "normal" Ly-break galaxies (Giavalisco, 
Steidel, & Macchetto 1996; Lowenthal et al. 1997). As 
might be expected from the previous section, the I-band 
images have even larger effective radii. All galaxies have 
a "red" nucleus, and the colors become bluer in the outer 
parts. 

Overall, the optical-to-infrared morphologies and sizes 
are strikingly similar to L* disk galaxies in the local uni- 
verse, with red bulges, more diffuse bluer exponential disks 
and scattered, UV-bright star forming regions. Some even 
show evidence of well-developed grand-design spiral struc- 
ture. However, the mean central surface brightness of the 
disks is 1 — 2 mag higher than that of nearby disk galaxies 
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and the mean rest-frame color (U — V) w is ~ 1 mag 
bluer (c.f. Lilly et al. 1998). Passive evolution can lead 
to disks with normal surface brightnesses at low redshift. 
Alternatively, the disks are disrupted later by interactions 
or evolve into SO's, which have higher surface brightnesses 
(e.g., Burstein 1979). 

5. SPECTRAL ENERGY DISTRIBUTION 

The overall SEDs of the galaxies show a large variety. 
Four of the galaxies (257, 267, 494 and 657) satisfy conven- 
tional U-dropout criteria (Madau et al. 1996; Giavalisco & 
Dickinson 2001). One galaxy is at too low redshift (302) 
to be classified as a U-dropout, and one other galaxy is 
too faint in the rest-frame UV (611). It has J s — K s > 2.3, 
and is part of the population of evolved galaxies identified 
by Franx et al (2003) and van Dokkum et al (2003). 

The red colors of the central components can be due 
to either dust, higher age, emission line contamination, 
or a combination of effects. The SEDs show stronger 
Balmer/4000 A breaks in the inner parts than in the outer 
parts (see Fig. 2 for an example). We derived colors in- 
side and outside of a 0.7 arcsec radius centered on the 
-fr s -band center. The mean differences are ~ 0.5 mag in 
observed 7si4 — Ks and ~ 0.2 mag in rest-frame U — V. 
Higher resolution NICMOS data are required to address 
the population differences in more detail. 
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Fig. 2. — The spectral energy distribution of source 611 in a 
Of! 7 circular diameter aperture (left) and in a concentric l" — 2" 
diameter ring (right), normalized to the Jf s -band flux. Over- 
plotted are independent model fits from Rudnick et al. (2003a). 

6. DISCUSSION 

We have found 6 large galaxies with characteristics sim- 
ilar to those of nearby disk galaxies: exponential profiles 
with large scale lengths, more regular and centrally con- 
centrated morphologies in the restframe optical than in 
the rest-frame UV, and, as a result, red nuclei. It is very 
tempting to classify these galaxies as disk galaxies, given 
the similarities with low redshift disk galaxies. However, 
kinematic studies are necessary to confirm that the ma- 
terial is in a rotating disk. Photometric studies of larger 
samples are needed to constrain the thickness of the disks. 
We note that simulations of Steinmetz & Navarro (2002) 
can show extended structures during a merging or accre- 
tion event. The expect duration of this phase is short, 
however, while the disk galaxies comprise a high fraction 
of the bright objects. It is therefore unlikely that a signifi- 
cant fraction of the galaxies presented here are undergoing 
such an event. 

The density of these large disk galaxies is fairly low: 
over a survey area of 4.7 arcmin 2 and to a magnitude 



limit of K s _tot = 22 they make up 6 out of 52 galaxies 
at 1.4 < z < 3.0. However, they do constitute 6 out of the 
12 most rest-frame luminous galaxies Ly > 6 x 1O 1o /iy O 2 L0 
in the same redshift range. The comoving volume density 
is ~ 3 x 10 hy Mpc~ 3 at a mean redshift z w 2.3. 
Obviously, larger area surveys are needed to establish the 
true density. We note that three of the galaxies only have 
photometric rcdshifts, one of which (a U-drop galaxy at 
Zphot — 1.82) is in the poorly tested range 1.4 < z < 2.0. 
The volume density of disk galaxies with r e > 3.6 h^} kpc 
in the local universe is much higher at <~ 3 x 10 h^ 
Mpc^ 3 (de Jong 1996), although many nearby disks would 
not be present in our high redshift sample because their 
surface brightness would be too low. 

We note that similar galaxies are absent in the very 
deep Near-IR imaging data on the HDF-N ( Williams et 
al. 1996; Dickinson 2000) Although notable differences be- 
tween optical and NIR morphologies were reported for two 
of the largest LBGs in the HDF-N, no large galaxies were 
reported to have red nuclei and exponential profiles as in 
the HDF-S. The two fields arc different in other aspects 
as well. We found earlier that the HDF-N is deficient in 
red sources (e.g., Labbe et al. 2003 versus Papovich, Dick- 
inson, & Ferguson 2001) and the disk galaxies are more 
luminous in the ii-band than most of the high-redshift 
galaxies found in the HDF-N (Papovich, Dickinson, & Fer- 
guson 2001). The larger number of red galaxies in the 
HDF-S (e.g., Labbe et al. 2003, Franx et al. 2003) and 
their strong clustering (Daddi et al. 2003) may indicate 
that the red galaxies and large disks are both part of the 
same structures with high overdensities, and evolve into 
the highest overdensities at low redshift, i.e. clusters. If 
this is the case, the large disk galaxies may be the pro- 
genitor of large SO galaxies in the nearby clusters, which 
have very similar colors as elliptical galaxies (e.g., Bower, 
Lucey, & Ellis 1992). 

Finally, we can compare the observed disk sizes to theo- 
retical predictions. It is often assumed (Fall & Efstathiou 
1980; Mo, Mao, & White 1998) that the disk scale length 
is determined by the spin parameter A and the circular ve- 
locity of the virialized dark matter halo (Fall & Efstathiou 
1980; Mo, Mao, & White 1998). For a ACDM cosmology, 
Mo, Mao, & White (1999) predict that the z ~ 3 space 
density of large (r e > 3.6 hj kpc) bright U-dropouts is 
1.1 xl0~ 4 h% Mpc- 3 , whereas we find ~ 2 x 10~ 4 h% 
Mpc~ 3 for our 4 U-drops at a mean < z >^ 2.4. This 
difference of a factor of two is not very significant, given 
our low number statistics and small survey volume. The 
combination of sizes and rotation velocities will give much 
stronger constraints on these models; it may be possible 
to measure the kinematics of some of these large galaxies 
using NIR spectrographs on large telescopes. 

We thank the staff at ESO for their hard work in tak- 
ing these data and making them available. This research 
was supported by grants from the Netherlands Foundation 
for Research (NWO), the Leids Kerkhovcn-Bosscha Fonds, 
and the Lorentz Center. GR thanks Frank van den Bosch 
and Jarle Brinchmann for useful discussions. 
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